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ABSTRACT 


The  evolution  of  a  two  dimensional,  incompressible,  rapidly  decelerating,  time- 
dependent  viscous  flow  about  a  sharp-edged  camber  is  simulated  through  the  use  of  the 
discrete  vortex  model.  Vorticity  is  represented  by  a  distribution  of  discrete  point 
vortices.  Each  vortex  is  convected  in  the  velocity  field,  calculated  locally  using  the  Biot- 
Savart  law.  The  roll-up  of  the  vortex  sheets,  the  distribution  of  velocity  and  pressure 
on  the  camber,  and  the  drag  force  are  calculated  at  suitable  time  intervals  for  a 
prescribed  time-dependent  flow.  Experiments  are  carried  out  in  a  vertical  water  tunnel 
partly  to  measure  the  drag  force  and  partly  to  record  on  a  video  tape  the  evolution  of 
the  v/ake.  The  measured  and  calculated  characteristics  of  the  flow,  such  as  the  growth 
of  the  wake  and  the  forces  acting  on  the  camber  are  found  to  be  in  good  agreement. 
Furthermore,  the  numerical  simulation  provided  a  plausible  explanation  for  the  cause 
of  parachute  collapse,  a  phenomenon  which  gave  impetus  to  the  numerical  and 
physical  experiments  described  herein.  The  numerical  model  developed  during  the 
course  of  the  investigation  is  applicable  to  any  time-dependent  flow  about  two- 
dimensional  cambered  plates  (circular  arcs  ). 
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I.  liNTRODUCTION 

A.       SEPAR.\TED  FLOWS 

The  separated  flow  about  blufT  bodies  has  been  almost  completely  unyielding  to 
both  analysis  and  numerical  simulation  for  a  number  of  mathematical  reasons  and 
fundamental  fluid  dynamic  phenomena.  Separation  gives  rise  to  the  formation  of  free 
shear  layers  which  roll  up  into  vortex  rings  or  counter-rotating  vortices.  They,  in  turn, 
interact  with  each  other,  with  the  counter-sign  vorticity  generated  at  the  base  of  the 
body,  and  with  the  motion  of  often  unknown  separation  points.  The  wake  becomes 
unsteady  even  for  a  steady  ambient  flow  and  the  problem  of  the  determination  of  the 
characteristics  of  the  wake  becomics  coupled  to  the  conditions  prevaihng  upstream  of 
the  separation  points.  Evidently,  viscosity  modifies  radically  the  inviscid  flow,  which, 
in  this  case,  cannot  serve  even  as  a  first  approximation  to  the  actual  flow.  The 
boundary  layer  equations  are  not  applicable  beyond  the  separation  points  and  are, 
therefore,  of  limited  use  in  bluff-body  flow  problems. 

The  separated  unsteady  flow  situations  involving  wake  return,  as  in  the  case  of  a 
decelerating  or  oscillating  body,  are  an  order  of  magnitude  more  complex  than  those 
where  the  vortices  continuously  move  away  from  the  body.  The  net  effect  of  the  wake 
return  is  twofold.  Firstly,  the  proximity  of  the  large  vortices  dramatically  affects  the 
boundary  layer,  outer  flow,  pressure  distribution,  and  the  generation  and  survival  rate 
of  the  new  vorticity.  Secondly,  the  vortices  not  only  give  rise  to  additional  separation 
points  and;  or  additional  vortices,  but  also  strongly  affect  the  motion  of  the  primary 
vortices.  These  effects  are  further  compounded  by  the  diffusion  and  decay  of  vortices 
and  by  the  three-dimensional  nature  of  the  flow. 

The  existing  finite  difference  and  finite  element  methods  cannot  yet  treat  the  high 
Reynolds  number  flows  with  sufficient  accuracy  for  a  number  of  reasons.  The  finite 
difference  schemes  require  a  very  fine  grid,  a  turbulence  model,  and  a  very  large 
computer  memory.  It  seems  that  the  modelling  of  the  turbulent  stresses  in  the  wake, 
particularly  in  time-dependent  flows  will  be  the  major  source  of  difficulty  in  all  future 
calculations.  Whether  or  not  it  will  ever  be  practical  to  apply  the  finite  difference  and 
finite  element  methods  to  high  Reynolds  number  flows  is  unknown.  The  inherent 
difficulties  are  certainly  significant  enough  to  warrant  exploring  other  solution 
methods. 
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Certain  separated  time-dependent  flows  may  be  simulated  through  the  use  of  the 
discrete  vortex  model  (DVM)  (see  e.g.,  Chorin  1973;  Sarpkaya  1975).  The  free  shear 
layers  which  emanate  from  the  sides  of  the  body  are  represented  by  an  assembly  of 
discrete  vortices.  The  strength  of  the  elemental  vortices  are  determined  through  the  use 
of  the  Kutta  condition.  The  use  of  a  suitable  convection  scheme  enables  one  to  march 
in  time  and  to  calculate  the  evolution  of  the  wake,  the  velocity  and  pressure 
distributions,  and  the  lift  and  drag  forces  acting  on  the  body.  The  work  described 
herein  deals  with  the  apphcation  of  the  DVM  to  decelerating  flow  about  a  two- 
dimensional  camber. 

B.       FLOW  ABOUT  A  CAMBERED  PLATE  AND  PARACHUTE  COLLAPSE 

The  determination  of  the  deployment  sequence  of  an  axisymmetric  porous 
parachute  and  the  unsteady  aerodynamic  loads  acting  on  it  present  a  very  complex 
coupled  problem.  The  development  of  an  analytical  or  numerical  model  which  takes 
into  account  the  effects  of  porosity,  gaps,  and  variable  opening  schemes  would  allow 
numerical  experiments  on  a  large  class  of  parachutes,  reduce  the  number  of  expensive 
field  tests  to  a  few  judiciously  selected  ones,  and  enable  the  designer  to  calculate  the 
time  history  of  the  fall  of  the  parachute  and  the  strength  required  to  survive  the 
aerodynamic  loads.  However,  the  development  of  such  a  model  is  hampered  by  a 
number  of  difficulties. 

The  previous  models  for  parachute  loads  are  based  by  and  large  on  empirical 
assumptions  (see  e.g.,  Heinrich  and  Saari  1987;  Mcwey  1972).  They  rely  on  the 
observation  that  families  of  parachutes  open  in  a  characteristic  length  and  seem  to 
have  aerodynamic  properties  that  relate  well  to  the  projected  area  of  the  parachute. 
The  apparent  mass  is  assumed  to  be  a  function  of  the  projected  area  only  and  is  not  a 
function  of  the  prevaihng  flow  characteristics.  The  vortex  sheet  analysis  was  used  by 
Klimas  (1977)  to  derive  the  acceleration-independent  apparent  mass  coefficient  for 
arbitrary-shaped  axisymmetric  surfaces.  Muramoto  and  Garrard  (1984)  used  a 
continuous-source  model  to  predict  the  steady-state  drag  of  ribbon  parachutes.  The 
analyses  did  not,  however,  deal  with  the  evolution  of  the  unsteady  wake  and  its 
interaction  with  the  canopy. 

It  is  in  view  of  the  foregoing  that  a  fundamental  study  of  the  separated  time- 
dependent  flow  about  two-dimensional  rigid  cambered  plates  was  undertaken.  Clearly, 
the  flow  about  a  rigid  cambered  plate  is  considerably  simpler  than  that  about  a  porous. 
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axisymmetric.  and  flexible  parachute  and  the  results,  regardless  of  the  degree  of  their 
agreement  with  corresponding  experiments,  may  not  have  direct  relevance  to  the 
practical  problem  under  consideration.  But  the  object  of  this  investigation  was  the 
understanding  of  the  evolution  of  the  wake  under  controlled  conditions  rather  than  to 
provide  a  design  tool.  It  is  hoped  that  an  investigation  of  this  type  will  reveal  the 
underlying  physics  of  the  phenomenon  (particularly  that  of  the  parachute  collapse), 
help  to  interpret  the  full-scale  results  and  will  provide  inspiration  for  the  development 
of  more  general  vortex  models  with  which  the  dynamics  of  axisymmetric,  porous,  and 
flexible  parachute  canopies  can  be  investigated. 
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II.  ANALYSIS 

A.       TR.ANS FORMATIONS  AND  THE  COMPLEX  VELOCITY  POTENTIAL 

The  calculation  of  the  velocity  of  any  one  of  the  vortices  and  the  force  acting  on 
the  body  requires  a  conformal  transformation  (in  which  the  camber  becomes  a  circle), 
a  complex-veiocity  potential  representing  the  vortices,  their  images,  and  the  two- 
dimensional  irrotational  flow  around  the  body,  and  the  use  of  the  generalized  Blasius 
theorem. 

The  flow  in  the  circle  plane  may  be  transformed  to  that  about  a  cambered  plate 
through  the  use  of  two  successive  transformations,  one  from  ^  plane  to  the  (^°  plane 
and  the  other  from  the  C,^  plane  to  the  z  plane.   These  are  given  by  (see  Fig.  2.1) 


(^-plane 


^•piane 


2ib 


-2ib 


Figure  2.1     Circle  and  physical  planes 
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z  =  ^«  -  —   ,       and  (,''=  (,  +  m  (2.1) 

Combining  the  two,  one  has  a  direct  transformation  from  the  C,  plane  to  the  z 
plane  as 

b2 

z  =  ^  +  m  -  (2.2) 

C,  +  m 

It  is  easy  to  show  that  the  camber  in  the  z  plane  is  a  circular  arc. 

The  y-axis  in  the  z  plane  passes  through  the  tips  of  the  camber.  It  is 
advantageous  to  locate  the  origin  of  the  coordinate  axes  at  the  geometric  center  of  the 
camber,  i.e..  at  the  center  of  the  circle  part  of  which  represents  the  camber.  This  is 
easily  accomplished  by  shifting  the  origin  of  the  coordinate  axes  by 

2nr  -  1 
m 

where  z  '  is  the  x  coordinate  of  the  orisin  of  the  circle  in  the  z  plane.   Thus,  one  has 


z  =C  +  m  -  -r  z    ,      with  z  '  =    -z  (2.4) 


which  transforms  the  circle  in  Fig.  2.1a  to  the  physical  plane  in  Fig.  2.  Id.    Table  1 
summarizes  the  relations! 
the  radius  of  the  camber. 


summarizes  the  relationship  between  m,  z  ',  the  included  angle  of  the  camber,  b,  and 


TABLE  1 
SUMMARY  OF  THE  PARAMETRIC  RELATIONSFIIP 

m                z  '          2a  b  R=  1/m 

o 

cos60  =  .5     -1     120  .866  2 

cos45  =  .707    0.     180  .707  V2 

cos30  =  .866   1/^3    240  .5  2/^3 
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The  complex   potential   function   W   in   tlie   circle   plane  (see   I"ig.    2.1a)   which 
describes  a  uniform  How  U  (assumed  to  be  time-dependent)  with  a  doublet  at  the 
origin  to  simulate  the  cylinder,  F^    clockwise-rotating  vortices  (called  q-vortices),  F. 
counter-clockwise  rotating  vortices  (called  p-vortices),  and  the  images  of  all  the  p-and 
q-vortices  in  the  circle  plane  may  be  written  as 

m    jp  m  jp  ^2  IP 

k.  =  r  k=  1  'kp 


ir  c^         '"    ir  "'    if  c^ 

2^  Qlq  k  =  l-^  '  k=r^  ^kq 


in  which  V.  and  ^j,  represent  respectively  the  strength  and  location  of  the  k-th  p- 
vortex,  Fj,  and  C,^  the  strength  and  location  of  the  k-th  q-vortex,  and  c  the  radius  of 
the  cylinder;  an  overbar  indicates  a  complex  conjugate.  The  need  for  the  separate 
identification  of  the  p-and  q-vortices  and  for  the  singling  out  of  one  of  the  vortices  in 
each  shear  layer  (namely  Fq  and  Fg  i.e.,  the  nascent  vortices  )  will  become  apparent 
later. 

B.       COMPLEX  VELOCITIES  OF  VORTICES 

The  convection  of  the  vortices  and  the  calculation  of  the  forces  acting  on  the 
body  require  the  evaluation  of  the  velocities  at  the  vortex  centers.  For  the  velocities  in 
the  circle  plane  this  reduces  to  subtracting  from  I:q.  (2.5)  the  complex  potential 
corresponding  to  the  vortex  for  which  the  velocity  components  are  to  be  determined 
and  evaluating  the  derivative  of  the  remaining  terms  at  C,  =  ^|^.  To  determine  the 
velocities  in  the  physical  plane,  however,  one  has  to  subtract  (iF, /2jr)  Ln(z  —  /,  )  from 
Eq.  (2.5)  or,  in  terms  of  (^  ,  the  terms  (see  e.g.,  Sarpkaya  1967,  1975) 


2k     ~'''^       ^""^        271  i;^ 


ir.  iF.  c^ 

^-Ln((;  -  r  )  +  — ^Ln(l   -  —  )  (2.6) 


k 
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It  should  be  noted  that  the  first  term  in  Eq.  (2.6)  is  the  complex  function 
corresponding  to  the  k-th  vortex  in  the  ^  plane.  The  second  term  appears  merely  as  a 
consequence  of  the  transformation  used.  It  is  easy  to  show  that  for  a  p- vortex  Eq. 
(2.6)  may  be  reduced  to  (see  e.g.,  Mostafa  1987) 

if,        (-b-)(L.     +  m) 
-  —JiP ^ -JM . (2  7) 

271      [{L,^^  +  m)-  +  b^  ]2 

A  similar  expression  can  be  written  for  the  q-vortex. 

C.       KUTTA  CONDITION 

The  fact  that  the  flow  separates  tangentially  with  a  finite  velocity  at  the  edges  of 
the  plate  (Kutta  condition)  may  be  expressed  by  requiring 


dW 

—    =  0  at   t;  =  i;    =-m  ±ib  (2.8) 


Thus,  inserting  Eq.  (2.5)  in  Eq.  (2.8)  one  has 

iFn.         1  1  ir.  1  1 


+  • 


JP  _  '  _      "  Oq 


ir,       1  1        ":Lir.        1 


2  i^  -._         >'         r  2 


k=i2^     ^t-^kp     r-jT      kti^JT      ^, -^kq     r-J. 
't    r  ^     r 

'kp  ^kq 


2 

-U(l-^)=0.  (2.9) 


Equation  (2.9)  may  be  decomposed  into  two  parts  as 
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'Op  ''Oq 

+  (-u^  +  ivj  =  0.  (2.10) 

where  the  terms  containing  the  strength  of  the  nascent  vortices  represent  the  velocity 
induced  at  the  tip  of  the  camber  by  the  nascent  vortices  and  the  term  in  parenthesis  the 
velocity  at  the  tip  due  to  all  other  vortices  (and  their  images),  the  doublet  at  the  center 
of  the  circle  in  the  C,  plane  and  the  ambient  velocity. 

Equation  (2.10)  represents  two  coupled  equations  for  the  strengths  and  positions 
of  the  nascent  vortices.  Thus,  the  solution  of  the  said  quantities  does,  in  general, 
require  an  iteration.  However,  this  iteration  may  be  avoided  by  noting  that  the 
velocity  induced  by  a  nascent  vortex  at  the  opposite  tip  is  very  small  and  certainly 
negligible.   Thus,  Eq.  (2.10)   for  one  of  the  nascent  vortices  may  be  reduced  to 


+  (-U  +iv  )  =  0  (2.11) 


^^oo            1          _ 

1 

2^          ^t"    ^Oq 

c2 

'^Oq 

A  similar  expression  may  be  written  for  the  other  nascent  vortex.  The  use  of  the  Kutta 
condition,  as  expressed  by  Eq.  (2.11),  will  be  further  explained  following  the  discussion 
of  the  tip  velocity.  It  suffices  to  note  that  all  nascent  vortices  satisfying  the  Kutta 
condition  do  not  yield  either  the  same  tip  velocity  or  the  same  velocity  distribution  in 
the  neighborhood  of  the  tip.  There  are,  in  fact,  certain  preferred  positions  for  the 
nascent  vortices  which  yield  physically  realistic  velocity  distributions  near  the  tips  of 
the  cambered  plate.  These  nascent  vortex  positions  have  been  determined  by  Mostafa 
(1987). 

D.       TIP  VELOCITY 

According  to  the  Kutta  condition  the  tangential  velocity  at  the  tip  is  finite.    It 
may  be  determined  through  the  use  of  I'Hopital's  rule. 
The  velocity  at  the  tip  is  given  by 


(2.12) 


dW     dW 

d; 

—  =  — 

at    z,  =  z    ±  2ib 

dz     di; 

dz 

t           o 
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The  use  of  THopital's  rule  then  yields, 


dW 
dz 


d^W     ib 

-2     {— )  (2.13) 


di;^        2 
z  =  z. 


Equation  (2.13)  yields  the  desired  finite  tip  velocities. 

E.       TIME  DEPENDENT-FORCES 

The  force  acting  on  the  body  in  the  physical  plane  may  be  calculated  either 
through  the  use  of  the  pressure  distribution  or  through  the  use  of  the  rate  of  change  of 
impulse. 

Bernoulli's  equation  for  unsteady  flow  is  given  by 


-> 


2 


p,        v.-  p,        V,-      j  dv 

(-^  +  -^  )  -  (  -^  +  — 1-)  -  j  —  ds  =  f  (t)  (2.14) 

p  2  p  2  J*'    ^t 


where  the  indices  indicate  two  points  on  the  body  in  the  physical  plane.  Since  there  is 
no  pressure  drop  across  the  shear  layer  and  since  the  integral  term  in  Eq.  (2.14)  is  zero 
at  the  tip  (i.e.,  ds=0),  one  has 

V   2        V   2 
f(t)  =  _^  -  _^  (2.15) 


where  V^^  and  V^  represent  the  tangental  velocities  on  the  upstream  and  downstream 
faces  of  the  tip.  It  is  important  to  note  that  f  (t)  in  Eq.  (2.15)  is  also  the  time  rate  of 
change  of  circulation,  i.e.,  the  rate  at  which  vorticity  is  shed  into  the  wake  from  the  tip 
of  the  cambered  plate. 

The  normalized  form  of  Bernoulli's  equation  between  any  two  points  m  and  n 
then  becomes 


p     _  p  V,--  V,2  V  -  -  V  -      5     "     V 

-^ — ^-^^  =  -^ ^-^  +  -2 :r-^  +—  f  2  T  ds  (2.16) 

pL  -2  U  ^  U  ^  ^t  •'       U  2 

•^     o  0  o  mo 
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The  integration  of  the  difTerential  pressure  between  the  upstream  and  downstream  faces 
of  the  camber  yields  the  force  components  in  the  x  and  y  directions,  i.e.,  the  drag  and 
lift  forces. 

The  force  acting  on  the  body  can  also  be  calculated  through  the  rate  of  change  of 
impulse.    It  is  given  by 

F  =  topc^U  (1  -  _j)  +  _[  r„  (z„  -  z„,)  1  (2.17) 

which  may  be  written  as 

F  cU         c  m^ 

C.  +  i  C,  = =  2k{ 5-)  { — )  (1  -  ^) 

^  >      2pU  -b         ^   U  2  ^     b  ^  2c- 

'^     o  o 


c      d         r.       C  c 

+  ^^n-  .>  .(T-^ta^)  -  f(^)]}  (2.18) 

2b  diht'c)    he         c  C 


in  which  U^  is  the  reference  velocity;  U,  the  rate  of  deceleration  of  flow  and  z=f^^j^), 
i.e.,  the  transformation  given  by  Eq.  (2.4).  Equation  (2.18)  may  also  be  deduced 
directly  from  the  generalized  Blasius  equation.  It  is  important  to  note  that  the  force 
calculated  from  Eq.  (2.18)  includes  the  effect  of  the  rate  of  change  of  circulation 
between  two  successive  time  steps.  Thus,  it  may  be  smaller  or  larger  (depending  on  the 
sign  of  F)  than  the  force  calculated  through  the  integration  of  the  instantaneous 
differential  pressure  Eq.  (2.16).  This  is  because  of  the  fact  that  the  instantaneous 
pressure  depends  only  on  the  prevailing  flow  conditions  and  does  not  account  for  the 
rate  of  change  of  total  circulation  between  successive  time  steps.  In  the  calculations  to 
follow  U    and  c  are  taken  as  unity  for  sake  of  simplicity. 

F.       METHOD  OF  CALCULATION 

The  methods  used  in  the  past  in  the  determination  of  the  vorticity  flux  from 
sharp-edged  bodies  may  be  roughly  classified  into  two  broad  categories.  The  first  of 
these  involves  the  use  of  variable  nascent  vortex  positions  (see  e.g.,  Sarpkaya  1968, 
1975)  and  the  second,  the  use  of  fixed  nascent  vortex  positions  (see  e.g.,  Clements 
1973-1975). 


In  the  present  analysis  the  method  of  variable  nascent  vortex  positions  is  used. 
To  explain  the  computational  details  of  the  method  let  us  consider  a  particular  time  t 
after  the  start  of  the  motion  and  assume  t  to  be  sufiiciently  large  so  that  there  are  a 
number  of  vortices  in  the  wake.  Then  the  appearance  and  convection  of  the  vortices 
proceed  as  follows: 

(1)  Select  a  vortex  position  along  the  radial  line  defined  by  6  =  117.72°  (see 
Mostafa  1987).    The  very  first  location  is  taken  r  =  1.1  . 

(2)  Calculate  the  strength  of  the  nascent  vortex  which  satisfies  the  Kutta 
condition.    This  is  an  exact  solution  and  requires  no  iteration. 

(3)  Place  the  nascent  vortex  at  the  corresponding  points  in  the  circle  and  physical 
planes  and  calculate  the  tip  velocity. 

(4)  Calculate  a  new  nascent  vortex  strength  from  0.5(V^^  -  V2")At  where  V^  is  the 
tip  velocity  and  V-,.  the  average  of  three  velocities  along  the  radial  line  in  the 
physical  plane,  i.e.,  at  r  =  1.05,  1.10  and  1.15. 

(5)  Compare  the  newly  calculated  circulation  with  that  obtained  from  the  Kutta 
condition.  If  the  diflerence  between  the  two  circulations  is  less  than  0.001 
procede  to  the  next  step.  If  the  said  difference  is  larger  than  0.001,  carry  out 
an  iteration  on  the  radial  location  of  the  nascent  vortex  as  many  times  as 
necessary  until  the  above  condition  is  satisfied.  If  the  circulation  calculated 
from  the  Kutta  condition  is  larger  than  that  calculated  from  the  tip  velocities, 
the  vortex  must  be  moved  towards  the  cylinder  and  vice  versa.  Also,  each 
time  the  direction  of  the  motion  of  the  nascent  vortex  is  changed  (inward  or 
outward),  the  marching  distance  is  halved  in  order  to  accelerate  the 
convergence  of  the  two  circulations. 

(6)  Calculate  the  velocity  induced  at  the  center  of  all  other  vortices; 

(7)  Convect  the  two  nascent  vortices  with  a  velocity  0.5(Vj  -I-  V-,  )  for  a  time 
interval  At  (note  that  the  vorticity  is  convected  with  the  average  velocity  of 
the  shear  layer). 

(8)  Convect  all  other  vortices  for  the  same  time-interval  At  using  a  second  order 
scheme  given  by 

z(t4-At)  =  z(t)  +  0.5  [3z(t)  -  z(t-At)  ]  At  (2.19) 

in  which  z  =  u  +  iv. 

(9)  Remove  the  vortices  from  the  calculation  whenever  they  come  nearer  than 
0.05  to  the  camber  in  the  physical  plane  (except  the  first  20  vortices  from  the 
tip); 

(10)      Coalesce  the  same  sign  vortices  with  a  separation  of  less  than  0.05  (in  the 
physical  plane,  except  the  first  20  vortices); 
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(11)  Calculate  the  tangential  velocities  and  pressures  on  the  inner  and  outer  faces 
of  the  camber.  Determine  the  drag  and  lift  forces  through  the  integration  of 
pressure  and  through  the  use  of  the  rate  of  change  of  impulse.  Make  plots  of 
suitable  variables  (e.g.,  velocity  distribution  near  the  tip,  variation  of  nascent 
vortex  circulation  with  time,  evolution  of  the  wake,  etc.); 

(12)  Check  the  flow  conditions  to  determine  the  state  of  the  calculations: 

(a)  If  V,  -  V,  >  0.1  repeat  the  foregoing  steps; 

(b)  Stop  the  introduction  of  nascent  vortices  if  0  <  Vj  -  V^  <  0.1  and 
return  to  step  Xo.  6; 

(c)  If  V2  >  Vj  switch  the  angular  positions  of  the  nascent  vortices  to  their 
image  points.  Calculate  V,  as  the  average  of  the  three  velocities,  at  the 
upstream  side  of  the  tip  of  the  camber,  at  three  radial  locations  (0.95, 
0.9,  and  0.85)  and  repeat  the  foregoing  steps;  and 

(13)  Make  plots  of  the  variations  of  various  flow  parameters  (e.g.,  tip  velocity, 
nascent  vortex  circulation,  evolution  of  the  wake,  force  coefficients,  etc.)  and 
terminate  the  run. 

The  foregoing  steps  are  quite  general  and  can  be  used  for  any  camber,  provided 

that  the  optimum  points  of  placement  of  the  nascent  vortices  are  determined  through 

the  use  of  a  method  developed  by  Mostafa  (1987). 
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III.  DISCUSSION  OF  RESULTS 

A.       NUMERICAL  AND  PHYSICAL  EXPERIMENTS 

The  calculations  were  carried  out  for  a  time-dependent  normalized  velocity  given 
by 


U  U    t 

—  =1      for    T  =  -^—    <  9.72  (3.1) 

L_  c 


and 


—  =  0.97T-0.05T2-3.70  (3.2) 

U 


and 

A  =  0.97-O.lOT  (3.3) 

in  the  interval  9.72  <  T  <  1 1.48  and 

U  =  -0.3423T    +  0.0072T2  +  3.82  (3.4) 

and 

A  =  -0.3423    +  0.01445T  (3.5) 

in  the  interval  11.48  <T.  <  17.95.  For  T  larger  than  18,  the  velocity  and  acceleration 
are  zero.  These  velocities  and  accelerations  correspond  to  that  encountered  in  a  series 
of  experiments  carried  out  in  a  vertical  water  tunnel.  A  detailed  description  of  the 
equipment  and  procedures  is  given  by  Sarpkaya  and  Ihrig  (1986)  and  will  not  be 
repeated  here.  Evidently,  the  calculations  can  be  carried  out  for  any  specified  variation 
of  the  velocity.    For  the  case  under  consideration,  the  fiow  begins  to  decelerate  at  T* 


=  9.72  and  the  velocity  of  the  ambient  flow  reduces  to  zero  at  about  T"  =  IS,  (see 
Fig.  3.1) 

The  computer  program  provided,  at  times  specified,  the  positions  of  all  the 
vortices,  the  rate  of  shedding  of  vorticity  fromi  the  tips  of  the  camber,  the  velocity 
distribution  on  the  upstream  and  downstream  faces  of  the  camber,  the  differential 
pressure  distribution,  and  the  force  coefficients. 

Figures  3.2  through  3.5  show,  at  T*  =  U  t/c  =  6.00,  the  evolution  of  the  wake, 
the  differential  pressure  distribution,  and  the  velocities  at  the  upstream  and 
downstream  faces  of  the  camber.  Figures  3.6  through  3.17  show  similar  plots  of  the 
wake,  pressure,  and  velocity  at  larger  times.  These  and  other  figures  (not  reproduced 
here  for  sake  of  brevity)  show  that  the  characteristics  of  the  flow  develop  symmetrically 
prior  to  the  onset  of  deceleration  (T*  <  9.72)  and  the  differential  pressure  is  positive 
everywhere  (i.e.,  the  pressure  inside  the  camber  is  larger  than  that  outside). 

Following  the  onset  of  deceleration  (see  e.g..  Figs.  3.8,  3.12,  and  3.16),  the 
differential  pressure  near  the  axis  of  the  camber  becomes  increasingly  negative.  The 
significance  of  this  result  is  that  had  the  model  been  flexible  (as  in  the  case  of  a 
parachute)  the  central  part  of  the  camber  would  have  collapsed  as  a  result  of  the 
particular  deceleration  it  is  subjected  to. 

For  T*  larger  than  about  13  (for  the  ambient  flow  under  consideration),  the 
velocities  induced  at  the  downstream  edges  of  the  camber  by  the  large  vortices  moving 
sideways  and  towards  the  camber  give  rise  to  oppositely-signed  vorticity.  This,  in  turn, 
leads  to  the  rapid  growth  of  the  secondary  vortices  (see  e.g..  Figs.  3.14  and  3.15  at  T" 
=  21.275).  Attention  is  drawn  to  the  fact  that  the  differential  pressure  shown  in  Fig. 
3.16  is  negative  over  a  large  central  portion  of  the  camber. 

Figure  3. IS  shows  the  variation  of  the  circulation  of  the  nascent  vortex  with  T*. 
As  expected,  the  vorticity  flux  is  quite  large  at  the  start  of  the  motion.  It  decreases 
gradually  prior  to  the  onset  of  deceleration  and  rapidly  thereafter.  For  T"  larger  than 
about  12,  i.e.,  when  the  secondary  vortices  come  into  existence,  the  circulation  becomes 
negative.  It  reaches  a  minimum  value  when  the  ambient  velocity  and  acceleration 
reduce  to  zero  (at  about  T*  =  18).  The  subsequent  motion  o[  the  primary  and 
secondar>'  vortices  increases  the  tip  velocity  and  hence  the  strength  of  the  nascent 
vortices.  Nevertheless,  the  circulation  remains  negative,  i.e.,  only  the  secondary' 
vortices  continue  to  receive  additional  vorticity. 
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Figure  3.19  shows  the  velocities  V^  and  V^  as  a  function  of  T*.  The  tip  velocity 
V^  decreases  from  an  initially  large  value  of  about  3  to  a  nearly  constant  value  of 
about  1.4  just  prior  to  the  onset  of  deceleration.  Subsequently,  V^  decreases  rapidly 
during  the  period  of  deceleration  and  prior  to  the  inception  of  the  secondar\' 
separation.  Then  V^  increases  to  about  2  because  of  the  backward  motion  of  the  large 
vortices  near  the  tips  of  the  camber.  Finally,  V^  decreases  once  again  as  the  primary' 
and  secondary  vortices  move  sideways  and  away  from  the  tips  of  the  camber  due  to 
their  mutual  induction  (see  Fig.  3.15). 

The  variation  of  V^  with  T*  is  significant  only  during  two,  relatively  short,  time 
intervals:  at  the  start  of  the  motion  and  at  the  start  of  the  deceleration.  These  are  the 
periods  during  which  the  vorticity  flux  changes  rapidly  in  order  to  maintain  the  Kutta 
condition.  During  the  remainder  of  time  V,  is  negligibly  small,  as  expected  on  the 
basis  of  the  pioneering  experiments  of  Fage  and  Johansen  (1928)  with  steady  flow  over 
various  types  of  bluff  bodies. 

Figure  3.20  shows  the  variation  of  the  drag  coelFicient  as  a  function  of  T".  It  is 
calculated  through  the  the  integration  of  pressure  about  the  camber.  Figure  3.20  also 
shows  that  C^  rises  rapidly  (due  to  the  rapid  accumulation  of  vorticity  in  the  growing 
vortices)  and  begins  to  decrease  as  the  vortices  develop  under  the  mfluence  of  a 
constant  ambient  velocity.  Then  the  force  decreases  sharply  at  the  onset  of 
deceleration  and  goes  through  zero  near  the  middle  of  the  deceleration  period  (T"  = 
14).  The  force  acquires  its  largest  negative  value  towards  the  end  of  the  deceleration 
period.    Subsequently,  the  force  gradually  decreases  to  zero. 

Figure  3.21  shows  a  comparison  of  the  calculated  and  measured  drag  coefllcients. 
In  seneral  the  agreement  between  the  calculated  and  the  measured  drag  coeilicient  is 
quite  good.  The  differences  are  primarily  due  to  the  fact  that  the  difTusion  of  vorticity 
has  not  been  taken  into  account  in  the  numerical  analysis.  It  is  possible  to  introduce  a 
small  artificial  reduction  in  circulation  in  order  to  bring  the  calculated  and  measured 
values  into  closer  agreement.  This  has  been  avoided  in  the  present  analysis  in  order  to 
keep  the  discrete  vortex  analysis  as  pure  and  simple  as  possible. 

Finally,  Fig.  3.22  shows  a  comparison  of  the  normalized  drag  forces  for  various 
periods  of  the  initial  steady  flow,  prior  to  the  onset  of  deceleration.  It  is  clear  that  the 
drag  forces  beyond  the  point  of  deceleration  are  nearly  identical.  In  other  words,  the 
force  acting  on  the  camber  is  not  materially  affected  by  the  duration  of  the  ambient 
steady  fiow  within  the  range  of  the  parameters  encountered  in  the  present  study. 
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Figure  3.1     Variations  of  the  velocity  and  acceleration  with  T* 
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Figure  3.2    Position  of  vortices  at  T*  =  6.00 
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Figure  3.3 


Instantaneous  velocity  of  the  vortices  at  T*  =  6.00 
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Figure  3.4    DilTerential  pressure  distribution  at  T*  =  6.00 
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Figure  3.5 


Velocity  distribution  on  both  faces  of  the  camber  at  T*  =  6.00 
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Figure  3.6    Position  of  vortices  at  T*  =  12.275 
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Figure  3.7    Instantaneous  velocity  of  the  vortices  at  T*  =  12.275 
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Figure  3.8 


Differential  pressure  distribution  at  T*  =  12.275 
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Figure  3.9    Velocity  distribution  on  both  faces  of  the  camber  at  T*  =  12.275 
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Figure  3.10     Position  of  vonices  at  T*  =  13.775 


o 
f-  ri 


O 
•d 


O 

m 

CO 


o 
■  oo 


••   .  '!• 


o 


o 

■  to 


o 


■  a 


X 


in       A 


o 


o 

I 


O 
■(JD 


O 
■  CO 


O 

-  o 


o 

tM 


0-8 


Q'9 


O'V 


0-^ 


O'O 

G/A 


O'S- 


0"r- 


0-9- 


o^a- 


35 


-\. 


/f/       --  \ 


^^f' 


%^^*^ 

'"^\ 


Hi 


*-^      » 


in 


a 
in 


o 


o 

■  en 


o 

'in 


o 

■  n 


d  > 


o 


I  I  I  1  I 

O^g       O'S       0">       Q"C       0'3 


OM 


I  i  I  I  I  I 

O'O       OM-     Q-S-     0-£-     0->-     0-3-  0"9- 


o 
■  ro 


o 

■  CO 

o 

CD  • 


0/A 


Figure  3.11     Instantaneous  velocity  of  the  vortices  at  T*  =  13.775 
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Figure  3. 


12    Differential  pressure  distribution  at  T*  =  13.775 
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Figure  3.13     Velocity  distribution  on  both  faces  of  the  camber  at  T*  =  13.775 
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Figure  3.14 


Position  of  vortices  at  T*  =  21.275 


a 
i-r\J 


o 
•d 


O 
in 


o 

"CO 


O 

i— 


o 

'U3 


o 

CM 


0-3 


0-9 


Q-f' 


O'S 


0/A 


Q-S- 


O't'- 


0-3- 


Q'O- 


Figure  3.15     Instantaneous  velocity  of  the  vortices  at  T*  =  21.275 
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Figure  3. 


16    Differential  pressure  distribution  at  T*  =  21.275 
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Figure  3.17    Velocity  distribution  on  both  faces  of  the  camber  at  T*  =  21.275 
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Figure  3.18    Variation  of  the  circulation  of  the  nascent  vortex  with  T* 
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Figure  3.19     Variation  of  the  velocities  V,  and  V^  with  T* 
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Figure  3.20    Calculated  drag  coefficient  as  a  function  of  T* 
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Figure  3.21     Comparison  of  measured  and  calculated  drag  coefficients 
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Figure  3.22    Measured  drag  cceriicienis  for  various  periods  of  initial  stead:/  flow 


47 


B.       CONCLUDING  REMARKS 

The  results  presented  in  this  section  have  shown  that  the  discrete  vortex  model 
can  be  used  with  confidence  to  predict  the  evolution  of  the  wake  about  a  cambered 
plate  immersed  in  an  arbitrary  time-dependent  flow.  The  drag  coefficients  resulting 
from  the  analysis  and  experiments  agree  reasonably  well.  This  agreement  can  be 
improved  with  the  introduction  of  a  small  circulation  dissipation. 

The  development  of  negative  difTerential  pressures  near  the  central  region  of  the 
camber  is  thought  to  be  primarily  responsible  for  the  inception  of  the  partial  collapse 
of  a  parachute  at  high  rates  of  deceleration.  This  phenomenon  takes  place  even  when 
the  total  drag  force  acting  on  the  parachute  is  still  positive.  The  sample  analysis 
presented  herein  also  shows  that  the  negative  differential  pressure  can  cover  a  large 
region  of  the  parachute  and  even  result  in  negative  drag.  The  basic  idea  emerging  from 
the  analysis  reported  herein  is  that  the  designs  which  incorporate  into  them  the  idea  of 
delaying  or  preventing  the  return  of  the  shed  vortices  to  the  canopy  (e.g., porosity 
management,  change  of  deceleration  history,  parachute  shape,  dissipation  and/or 
destruction  of  the  organized  wake)  will  be  the  ones  which  could  avoid  the  collapse 
phenomenon.  Extensive  analysis  and  small  scale  experiments  coupled  with  a  few 
judiciously  selected  field  tests  may  help  to  arrive  at  practically  and  phenomenologically 
sound  parachute  designs. 
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